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Abstract 

The problem of finding, in an edge- weighted bidirected graph G = [V,E), a cycle 
with minimum mean weight of its edges generalizes similar problems for both directed 
and undirected graphs. (The problem is considered in two variants: for the cycles 
without repeated edges and for the cycles without repeated nodes.) In this note we 
develop an algorithm to solve this problem in 0{V'^ min(F^, Slog V^))-time (to compare: 
the complexity of an improved version of Barahona's algorithm for undirected cycles 
is 0(F^)). Our algorithm is based on a certain general approach to minimum mean 
problems and uses, as a subroutine, Gabow's algorithm for the minimum weight 2-factor 
problem in a graph. The problem admits a reformulation in terms of regular cycles in 
a skew-symmetric graph. 
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1 Introduction 

Among a variety of discrete optimization problems, an interesting class is formed by the 
problems consisting in finding an object (subset) of a given type in which the mean weight 
of an element is minimized. One of the most popular problems of this sort is the minimum 
mean cycle problem in an edge-weighted directed graph. The classical algorithm due to 
Karp jni, based on a dynamic programming approach, finds such a cycle in 0{nm) time. 
(Here and later on n and m denote the numbers of nodes and edges, respectively, in the 
input graph.) An analogous problem for cycles of an undirected graph can also be solved 
efficiently: Barahona reduces it to a series of minimum weight 0-join computations (see 
also (HI Sec 29.11]), and an improved version of his algorithm requires only 0{n) such 
computations. 
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This note presents an algorithm with complexity 0(n^min(n^,mlogn)) for a more gen- 
eral problem, namely, for finding a minimum mean cycle in an edge-weighted bidirected 
graph. (The concept of bidirected graph was introduced by Edmonds and Johnson 0; more 
about bidirected graphs can be found in, e.g., jH|.) 

Recall that in a bidirected graph G = {V, E) edges of three types are allowed: a usual 
directed edge, or an arc, that leaves one node and enters another one; an edge directed 
from both of its ends; and an edge directed to both of its ends. When both ends of an edge 
coincide, the edge becomes a loop; to slightly simplify our description we admit no loop 
entering and leaving its end node simultaneously (this will lead to no loss of generality in 
what follows). 

A walk in G is an alternating sequence P = (s = vq, ei, wi, . . . , Cfc, ffc = t) of nodes 
and edges such that each edge Cj connects nodes Vi-i and fj, and for i = 1, . . . , /c — 1, the 
edges Cj, Cj+i form a transit pair at u,, which means that one of Cj, Cj+i enters and the other 
leaves Vi. (Note that ei may enter s and may leave t.) If, in addition, k > \, vq = Vk 
and the pair ei,efc is transit at uq, P is called a cycle. A cycle is called node- simple if all 
its nodes (except vq = Vk) are different, and edge-simple if all its edges are different. Note 
that G may contain edge-simple cycles but no node-simple ones. 

Our main problem is stated as follows: 

(P) Given a bidirected graph G = (V, E) and a function w : E —> M of weights of edges, 
Find an edge-simple cycle G of G that minimizes the value w{G) := w{G)/\G\. 

Here w{C) denotes the total weight of edges, and \G\ the number of edges (the length) of 
C. So we deal with the minimum mean edge-simple cycle problem in an edge-weighted 
bidirected graph. Instead, one can state, in a similar way, the minimum mean node-simple 
problem in (G, w), which also may find reasonable applications. However, the latter is easily 
reduced to (P) (and therefore, can be excluded from consideration in what follows). Indeed, 
modify the graph as follows: replace each node f by a pair of nodes ui, f2, connecting them 
by an edge with zero weight from vi to V2, and each edge of the original graph that enters 
(leaves) v make entering vi (resp. leaving ^2). The obtained graph has 2n nodes and m + n 
edges, all edge-simple cycles G' in it are node-simple and just correspond to the node-simple 
cycles G of the original graph; moreover, for corresponding G' and C, one has w{G') = w{G) 
and \G'\ = 2\G\. 

In fact, the minimum mean cycle problem for directed graphs G becomes equivalent to 
(P) (as well as to its node-simple version) when G is considered as a bidirected graph. As for 
the undirected case, an undirected graph G can be turned into a bidirected one by assigning 
the orientation of each edge from both of its ends and by adding, for each node v, a loop 
with zero weight that enters v (twice). This reduces the undirected minimum mean cycle 
problem to problem (P) (but not to its node-simple version). 

We show the following. 

Theorem 1.1 Problem (P) can be solved in 0(n^ min(n^, m log n)) time. 

Another class of nonstandard graphs is formed by so-called skew-symmetric graphs — 
directed graphs with involutions on the nodes and on the arcs that change the orientation 
of each arc. (This class of graphs, under the name of antisymmetric al digraphs, was intro- 
duced by Tutte jH].) There is a close relationship between these and bidirected graphs, and 
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problem (P) is, in fact, equivalent to the problem of finding a regular cycle having minimum 
mean weight in a skew-symmetric graph with symmetric weights of arcs (precise definitions 
will be given later). As a consequence, we obtain an 0(n^ min(n^, m log n))-algorithm for 
the latter problem. 

A natural approach to the minimum mean problem for an abstract family T of subsets 
of a set consists in reducing this problem to a series of (auxiliary) problems of finding a 
member of T with minimum total weight when the weight function is shifted by a constant. 
The shifts are chosen in such a way that the cardinalities of intermediate subsets in this 
process monotonically decrease. Therefore, the process is finite, provided that the auxiliary 
problem is well-solvable. In a more general approach, T is extended to a family V by adding 
certain sets represented as the disjoint union of members of T . A clever choice of T> may 
simplify the auxiliary problem significally, yielding an efficient algorithm for the original 
problem. Just this idea is applied in ^ where undirected cycles are extended to 0-joins. 
(One more approach to a certain class of minimum mean problems and related topics are 
discussed in 0.) 

In our case, the role of extended family T) plays the family of circulations that traverse 
every node of a given bidirected graph at most twice. Then the auxiliary problem becomes 
equivalent to the minimum weight 2-factor problem in a certain associated undirected graph. 
The latter can be solved in 0(nmin(n^,mlogn)) time by use of Gabow's algorithm 1^. This 
yields an algorithm with the desired complexity, taking into account that the number of 
iterations in the process is 0(n). 

The above-mentioned general approach is described in Section [21 and a proof of Theo- 
rem ^3 is given in Sectional The related problem for skew-symmetric graphs is discussed 
in Sectional 

2 A General Approach 

In a general setting, the minimum mean problem is formulated as follows: 

(M) Given a family T of nonempty subsets of a finite set E and a function u; : ^ R, 
Find X ^ T minimizing w{X) := 'w{X)/\X\, 

where w{X) denotes X^eex ^(fi)- For convenience the elements of ^ are called feasible sets; 
so the goal is to find a feasible set with minimum mean weight. 

Let us call a family P of nonempty subsets of -E a disjoint extension oi J- \i J- C^T) and 
each member X oiV can be represented as the disjoint union of feasible sets, i.e., 

(1) X = Xi U . . . U Xk for pairwise disjoint sets Xi, . . . , X^ G J^. 

(To emphasize that the sets involved in the union are pairwise disjoint, we use notation U 
rather than U.) Let s(V) denote the maximum cardinality (size) of a member of T>. 

For an arbitrary disjoint extension V of problem (M) can be reduced to at most s(I?) 
problems of finding a set X G P having minimum weight w'{X) with respect to some other 
weight function w' , followed by one problem of constructing a decomposition Q for the 
final set X. The idea is rather transparent and has been encountered in special cases (cf., 

e.g., 
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Functions w' occurring in the desired reduction are obtained by shifting it; by a constant. 
More precisely, we say that w' is the shift of w by a G M if w'{e) := w{e) — a for all e € E. 
We also denote w' by w"' and call a the shift number. The method below relies on two easy 
properties: 

(2) w'^{X) = w{X) — a for any nonempty subset X C E; 

(3) if X = Xi U . . . U Xk for nonempty subsets Xj, then w{X) > miuj uJ(X'j). 

Method. In the beginning put a := max(u;(e): e S E). Then iteratively construct a set X 
and update a as follows. At each iteration, take as X a member of V with w"'{X) minimum, 
by solving the corresponding minimum weight problem. Put b := w'^{X). If (i) 6 = 0, then 
the process terminates with constructing a decomposition Q of X and outputting any of 
its members. Otherwise (ii) update a := a + b and proceed with the next iteration. 

Lemma 2.1 The above method terminates in at most s{V) + 1 iterations and outputs an 
optimal solution to problem (M). 

Proof. 

First of all we observe that at each iteration there exists a set y G P satisfying w"'{Y) < 
0. Indeed, the initial choice of a guarantees this property to hold at the first iteration. And 
for each iteration, we have w"" (X) = w''{X) - 6|X| = 0, where a is the shift number at the 
beginning of this iteration, X is the set in V found on it, 6 = IZJ'*(X), and a' is the next shift 
number a + b. Then the property holds for the next iteration. 

Thus, each number b in the process is nonpositive. Suppose 6 = happens at some 
iteration. Then the set X found on this iteration satisfies w^{X) = 0, and by the minimality 
of X, we have w^iY) > for all Y e V; equivalently: = ^''(X) < w^iY). Now and 
imply that W{Xi) = a for each member X, in a decomposition of X, and that W(Y) > a for 
any Y G (taking into account that .7-" C P). So Xj is an optimal solution to (M). 

Finally, consider two consecutive iterations. Let a and X denote the shift numbers at 
the beginning of the former iteration and the set found on it, and let a' and X' denote similar 
objects on the latter one. Suppose that the latter iteration was not the last in the process; 
then W'iX') < 0. Also a' = a + uJ'^(X) imphes W' {X') = uJ"(X') - tU«(X). Therefore, 

w;"(X')/lX'[ < w^X)/\X\. 

This inequality is possible only if |X'[ < |X|, in view of w"(X) < w'^{X') and w"-{X) < 0. 
Thus, the process is finite and the number of iterations does not exceed s{'D) + 1. □ 

Remark. It is seen from the above proof that the number of iterations in the method 
is estimated via the minimum cardinality of a set X £ 7? minimizing w°'{X) on the first 
iteration. Also the method can start with any initial shift number a for which one guarantees 
the existence of a set y € P with w"'{Y) < (provided that such an a can be computed 
efficiently) . 

To illustrate the method, consider problem (M) for the family of simple cycles, or 
circuits, in a undirected graph G = (V, E) with edge weights w (regarding circuits as edge 
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sets). Since the minimum weight problem for circuits is NP-hard when negative weights are 
possible, and therefore, P := ^ is a bad choice, one has to take as V some nontrivial disjoint 
extension of T. Following jT|, we put V to be the family of nonempty 0-joins in G. (Recall 
that for T C y with \T\ even, a set J C is called a T-join if the set of nodes with odd 
degrees in {V, J) is exactly T.) Clearly a 0-join is decomposed into pairwise (edge-) disjoint 
circuits; such a decomposition is constructed in linear time. Each iteration of the method 
in our case consists in finding a minimum weight 0-join, which can be carried out in O(n^) 
time (see 0). Also the number of iterations is 0{m). Hence the obtained algorithm runs 
in 0{n^m) time. 

In an improved version of this algorithm, the number of iterations reduces to 0{n), 
yielding the overall time bound 0{n^). This relies on the existence (and the possibility to 
efficiently construct) an initial shift number a such that there exists a nonempty 0-join J 
whose weight w'^{J) is minimum and whose cardinality is 0{n), and moreover, this weight 
is nonpositive; cf. Remark above. 

To compute the desired a, we order the edges of G by nondecreasing their weights: 

w{ei) < w{e2) < ■■■ < w{em)- 

We find the minimum number k such that the graph Gk = {V, {ei, . . . , ek}) contains a 
circuit C. Put a := w{C); then w'^iC) = 0. Define 

:= {eeE\ w"(e) > 0} 
E- := {eeE \ w"(e) < 0} 

Let J be a nonempty 0-join with \J\ minimum among all 0-joins J' minimizing w°'{J'). The 
circuit C is also a 0-join, hence w'^{J) < w°'{C) = 0. We assert that |J| < 2n. Indeed, 
since Gfc-i is acyclic and, obviously, contains E' ^ we have \ J r\E^\ < \E~\ <k <n. Also 
I J n E^\ < n. For otherwise the subgraph {V, J PI E^) would contain a circuit C . Then 
J' := J — C is a nonempty 0-join satisfying w°'{J') < w"'{J) and \ J'\ < [ J[, contradicting 
the minimality of J. Thus, | J| < 2n, and a is as required. 

3 Minimum Mean Cycles in Bidirected Graphs 

In this section we specialize the method described in Section |21 to find a minimum mean 
edge-simple cycle in a bidirected graph G = {V, E) with a weighting w: E ^ on the 
edges. For brevity we omit the adjective "edge-simple" in what follows. 

Let us call a subset X Q E balanced if for each node v, the numbers of edges in X 
entering v and leaving v are equal (counting twice each doubly entering or doubly leaving 
loop at V, if any). In particular, the edge-set of any cycle is balanced. If each node is entered 
by at most two edges in X, we say that this balanced set is small. In a similar fashion, a 
cycle is called small if it passes each node at most twice, i.e., its edge-set is a small balanced 
set. Balanced sets and small cycles are related by the following property (in fact, known 
in literature), which is a bidirected analog of the fact that a circulation in a digraph is 
decomposed into simple cycles. 

Lemma 3.1 Each balanced set X is representable as the union of pairwise edge-disjoint 
small cycles (regarding a cycle as an edge set). 
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Proof. 

Assuming X ^ %, choose a maximal edge-simple walk P = (t>o, ei, fi, . . . , ek,Vk) in the 
subgraph of G induced by X. The fact that X is balanced easily implies that vq = Vk and 
that ei, Cfc form a transit pair at uq, i.e., P is a cycle. 

Deleting the edges of P from X, we again obtain a balanced set. So it follows by 
induction that X is representable as the union of pairwise edge-disjoint cycles. Let the 
number of cycles in such a representation be as large as possible. We assert that each cycle 
is small. 

Indeed, suppose some cycle C = {vo,ei,vi, . . . ,ek,Vk) among these is not small, i.e., C 
passes some node v at least three times. Since a cycle can be considered up to reversing 
and cyclically shifting, we may assume that v = vq = Vi = Vj for some < i < j < k, and 
that ei leaves vq. If for some p G {hj}, the edge Cp enters Vp, then the pair ei,ep (as well 
as efc,ep+i) is transit at v, and we can split C into two cycles: the part of C from vq to Vp 
and the rest. And if Cj leaves Vi and ej leaves Vj, then the pair Cj+i, ej is transit, and C can 
again be split into two cycles. A contradiction. □ 

It is straightforward to devise a linear time algorithm that decomposes a balanced set 
into small cycles. 

An immediate consequence of ^ and Lemma IXTl is that for any cycle C, there exists 
a small cycle C such that w{C') < w{C). Therefore, instead of all cycles, we can take 
as the family T in problem (M) the collection of (the edge-sets of) small cycles. Also the 
desired disjoint extension V oi can be defined to be the collection of small balanced sets. 
The cardinality of any small balanced set X does not exceed 2n since each node is incident 
with at most four edges in X. This implies that the number of iterations of the method 
applied to these T and V is 0{n), by Lemma ITTl Each iteration consists in finding a small 
balanced set X with minimum weight w'{X) for a current weight function : ii^ — > M, and 
it remains to explain how to solve this problem. 

We reduce it to the minimum weight 2-factor problem in a certain undirected multigraph 
G = {V,E), with possible loops, formed from G as follows. Each node v €V generates two 
nodes vi,V2 in G. Each edge e in E connecting nodes u and v generates an edge e with the 
same weight connecting nodes Ui and Vj, by the following rule: if e enters u then i = 1, 
otherwise i = 2, and similarly for v and j. In particular, a doubly entering loop at v (if 
any) induces an undirected loop at vi. Finally, for each v & V, we add two parallel edges 
connecting vi and V2 in G and assign them zero weight; these edges are called auxiliary. An 
example of this transformation is depicted in Fig. ^ 

Recall that a 2-factor in a (multi)graph is a subset of edges such that each node of the 
graph is covered by exactly two edges of this subset (counting a loop twice). There is a 
natural correspondence cj) between the small balanced sets in G and the 2-factors in G (up 
to swapping parallel auxiliary edges), and this correspondence preserves set weights. More 
precisely, for a small balanced set X in G, (p{X) contains all edges of G generated by X. In 
addition, for each node v £ V incident with d entering edges in X (counting a loop twice), 
(j){X) contains 2 — d auxiliary edges connecting vi and V2- One can check that (j){X) is 
indeed a 2-factor. Moreover, any 2-factor of G can be obtained in this way. 

Applying the 0(n min(n^, m log n))-algorithm due to Gabow j2] to find a minimum 
weight 2-factor in G, we conclude that each iteration in our method can be performed 
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Figure 1: Reduction to 2-factor problem. 

with a similar time bound. Since the number of iterations is 0(n), we obtain an algorithm 
to solve the original problem (P) in 0(n^min(n^,mlogn)), as required in Theorem ll.il 

4 Skew- Symmetric Graphs 

The above result on bidirected graphs can be translated into the language of skew-symmetric 
graphs. (About skew-symmetric graphs and problems on them, see, e.g., (HIEIEI-) 

A skew- symmetric graph is a digraph G = {V,E), with possible multiple arcs, endowed 
with two bijections av,crE such that: ay is an involution on the nodes (i.e., crv{v) ^ v and 
av{(Tv{v)) = V for each node v); cfe is an involution on the arcs; and for each arc a from u 
to f, cT£;(a) is an arc from cry(w) to av{u). For brevity, the mappings av,crE are combined 
into one mapping a on ^ U ii^, which is called the symmetry (rather than skew-symmetry) 
of G. For a node (arc) x, its symmetric node (arc) a{x) is also called the mate of and we 
will use notation with primes for mates, denoting a{x) by x'. 

Observe that if G contains an arc a from a node v to its mate v' , then a' is also an arc 
from V to v'. 

The symmetry a is extended in a natural way to walks, cycles and other objects in 
G. In particular, two walks or cycles are symmetric to each other if the elements of one 
of them are symmetric to those of the other and go in the reverse order: for a walk P = 
{vQ,ai,vi,. . . ,ak,Vk), the symmetric walk a{P) is (w^, a'^, . . . ,a[,v'Q). 

Next we explain a relationship between skew-symmetric and bidirected graphs (cf. 
Sec. 2]). Given a skew-symmetric graph G = {V,E), choose an arbitrary partition vr = 
{^1,^2} of V such that V2 is symmetric to Vi. Then G and vr determine the bidirected 
graph G = {V, E) with V := V\ whose edges correspond to the pairs of symmetric arcs in 
G. More precisely, arc mates a, a' of G generate one edge e of G connecting nodes tt, t> E Vi 
such that: (i) e goes from u to f if one of a, a! goes from utov (and the other goes from v' 
to m' in V2); (ii) e leaves both u^v if one of a, a' goes from u to v' (and the other from v to 
w'); (iii) e enters both u^v if one of a, a' goes from u' to v (and the other from v' to u). In 
particular, e is a loop if a, a! connect a pair of symmetric nodes. 

Conversely, a bidirected graph G = iV , E) determines a skew-symmetric graph G = 
iy, E) with symmetry a as follows. Take a copy a{v) of each element v of V, forming the 
set V := {a{v) \ v G F|. Now set V := V U V . For each edge e of G connecting nodes u 
and f , assign two "symmetric" arcs a, a' in G so as to satisfy (i)-(iii) above (where u' = a{u) 
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(a) Bidirected graph G. 



(b) Corresponding skew- 
symmetric graph G. 



Figure 2: Related bidirected and skew-symmetric graphs. 

and v' = (j{v)). An example is depicted in Fig. |21 

(Note that one bidirected graph generates one skew-symmetric graph, by the second 
construction. On the other hand, one skew-symmetric graph generates a set of bidirected 
graphs, depending on the partition vr of by the first construction. One can check that 
among these bidirected graphs, one is obtained from another by the following operation: 
choose a subset X of nodes, and for each node v m. X and each edge e incident with v, 
reverse the direction of e at v. This implies that the sets of walks (cycles) in these bidirected 
graphs are the same.) 

There is a one-to-one correspondence between walks in G and (directed) walks in G. 
More precisely, let r be the natural mapping olV U E to V U E (obtained by identifying 
the pairs of symmetric nodes and arcs). Each walk P = (wq, ai, wi, • • • , flfc, ^fc) in G (where 
tti is an arc from vi-i to Vi) induces the sequence 

t{P) ■= (T(wo),r(ai),r(t!i),...,r(afc),r(t)fc)) 

of nodes and edges in G. One can see that t{P) is a walk in G and that t(P') is the walk 
reverse to r(P). Moreover, for any walk P in G, there is exactly one preimage t"^(P). 

Finally, following terminology of j^, a walk (cycle) in a skew-symmetric graph is called 
regular if it is arc-simple and contains no pair of symmetric arcs (while pairs of symmetric 
nodes in it are allowed). Certain problems on regular walks and cycles are studied in OE] 
(and some other works). One more problem is: 

(S) Given a skew- symmetric graph G = (V, E) and a symmetric weight function w : E 

(i.e., w{a) = w{a') for all a & E), find a regular cycle C in G whose mean weight 
w{G) is minimum. 

It is not difficult to check that for each regular cycle C in G, the corresponding cycle 
r(C) in the bidirected graph G is edge-simple, and vice versa. Moreover, the weights (and 
the mean weights) of the corresponding cycles are equal, assuming that for e € E, the weight 
of r(e) is defined to be ?/;(e) (this is well-defined since w is symmetric). Therefore, problem 
(S) is equivalent to (P), and we obtain the following. 
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Corollary 4.1 Problem (S) can he solved in 0(n^ min(n^, m log n)) time. 
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